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ABSTRACT 

We report the first detection of a redshift-depth enhancement of background galaxies magnified by 
foreground clusters. Using 300 000 BOSS-Survey galaxies with accurate spectroscopic redshifts, we 
measure their mean redshift depth behind four large samples of optically selected clusters from the 
SDSS surveys, totalling 5 000 — 15 000 clusters. A clear trend of increasing mean redshift towards 
the cluster centers is found, averaged over each of the four cluster samples. In addition we find 
similar but noisier behaviour for an independent X-ray sample of 158 clusters lying in the foreground 
of the current BOSS sky area. By adopting the mass-richness relationships appropriate for each 
survey we compare our results with theoretical predictions for each of the four SDSS cluster catalogs. 
The radial form of this redshift enhancement is well fitted by a richness-to-mass weighted composite 
Navarro-Frenk- White profile with an effective mass ranging between M200 ~ 1.4-1.8 x 10 14 M Q for 
the optically detected cluster samples, and M 2 oo ~ 5.0 x 10 14 M Q for the X-ray sample. This lensing 
detection helps to establish the credibility of these SDSS cluster surveys, and provides a normalization 
for their respective mass-richness relations. In the context of the upcoming bigBOSS, Subaru-PFS, 
and EUCLID-NISP spectroscopic surveys, this method represents an independent means of deriving 
the masses of cluster samples for examining the cosmological evolution, and provides a relatively 
clean consistency check of weak-lensing measurements, free from the systematic limitations of shear 
calibration. 

Keywords: galaxies: clusters: general — gravitational lensing: weak — cosmology: observations — 
dark matter 



1. INTRODUCTION 

Our current paradigm predicts that the properties of 
dark energy will be imprinted in the matter power spec- 
trum and its time evolution. Although probing the large 
scale structure in the universe is a promising strategy to 
investigate the different dark-energy model candidates, 
it is becoming increasingly challenging to detect the rel- 
atively subtle signatures of the various lesser ingredients 
and higher order effects, forcing us to observe many in- 
dependent physical phenomena over gigantic volumes, to 
reach sufficient joint statistical precision. 

In the era of large scale cosmological surveys, gravi- 
tational lensing is one of the most promising methods 
for establishing the fluctuations of matter in the uni- 
verse, as lensing is purely gravitational and can be largely 
assumed to be independent of the unknown nature of 

the dark matte r. One may employ two-dimensional 

(|Kilbinger et al.ll2013|) or tomographic ()Beniamin et al.l 
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|20~L1 cosmic shear to measure the integrated lensing ef- 
fect of matter flu ctuations along the line-of-sight up t o 

[2001 
2010) 



very large scales ([Taylor et al.11200-3: IMassev et al.l 



Or, one may use cluster abundances (|Rozo et al 



to put constraints on the evolution of the dark matter 
halo mass function, which is known to be a very sensitive 
function of the evolution of the mean mass density and of 
the dark energy properties. Masses can be derived from 
the radial s hear profile around c l usters in the weak lens- 
ing regime (ISheldon et all 120091; ILeauthaud et al.l 120101: 
IQguri fe HamanallgOllHUmetsu et al.ll2012fl. from X -ray 



(jRvkoff et al.ll2008| ). SZ effect (iMarriaee et al.ll2011| ). or 
satellite kinematics (|More et al.ll201lf ). 

Even though lensing from stacked clusters undoubtedly 
leads to a less noisy mass estimate among all these meth- 
ods, accurate statistics will be limited, however, by our 
ability to measure galaxy shapes down t o the required 
precision. Recently. iHevmans et al.l <|2012f ) demonstrated 
that for a survey like the CFHTLS-wide, spreading over 
154 deg 2 , the level of systematics in shape measurement, 
as given by the state-of-the art sha pe measurement tech- 
nique (Lensfit, iMiller et alJ[2uTl . was lower than sta- 
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tistical errors. For upcomin g experiments such a s the 
Hyper Suprime Cam fHSC. iMivazaki et all I2012D sur- 
vey, the Dark Energy Survey (DE33), the Large Synoptic 
Survey Telescope pr oject (LSST0) or EUCLIE0 survey 
(jLaureijs et al.11201 lh which will probe volumes orders of 
magnitude larger than the CFHTLS-wide, the limiting 
requirements will be much more stringent. 

Lens magnification provides independent observa- 
tional alternatives and complementary mea ns to 

weak shape measurem ent (jVan Waerbeke et al.1 120101 : 
iHildebrandt et al~ll2011h . In practice magnification bias 
has been used to improve the mass profiles of individual 
massive clusters, and stacked samples of clusters from re- 
cent dedicated cluster surveys. Here the number counts 
of flux limited samples of background sources is reduced 
in surface density by the magnified sky area, and com- 
pensated to some ex tent by fainter objects magnified 
above the flux limit (|Broadhurst et aD Il995l hereafter 
BTP95) with a net magnification-bias in the surface den- 
sity depending on the slope of the background counts. 
For very red background galaxies behind clusters the net 
effect is a clear "depletion" of red background galax- 
ies to wards the cluster center (jBroadhurst et al.l I2005L 
2008) and is used in combination with weak shear to 
enha nce the precision of individual cluster mass pro - 
files (jUmetsu fc Broadhursti[2008l; lUmetsu et al.ll2011bD . 
A strength of this method resides in the fact that no 
shape information is used and is therefore less prone 
to systematic errors, as opposed to shear measurement 
for which complex corrections are required (jKaiser et al.l 
[1991 lUmetsu fc Broadhurs^[200l 

In the statistical regime, for deep wide area surveys, 
progress has been made recently in a similar way us- 
ing distant Lyman-break galaxies behind stacked sam- 
ples of foreground clusters, resulting in a claimed pos- 
itive magnification-bias enhancement of the observable 
surface density of background galaxies at small radius. 
In practice, some cross contamination of cluster mem- 
bers and foreground galaxies with the background will 
complicate such me asurements. Unti l recently, magnifi- 
catio n bias studies (|Ford et al.l T2012: H ildebrandt et al.l 
|2013[ ) have faced such a high statistical error that their 
conclusions were not affected by these systematics. How- 
ever in the future a careful treatment of potential sources 
of systematics will be mandatory. 

With redshift information for large samples of back- 
ground galaxies, cluster magnification may be tackled 
more fully, as the luminosity function of background 
galaxies is magnified both in the density and luminos- 
ity directions, in a characteristic redshift dependent way. 
For flux limited redshift samples, the redshift distribu- 
tion becomes skewed to higher mean redshift as a conse- 
quence of magnification, and the mean luminosity is also 
modified owing to the curvature (non-power law) of the 
luminosity function (BTP95). The latter effect has been 
claimed in the case of magnified QSO's behind galaxies 
in the SDSS survey (jMenard et alJhoiOft . 

In this paper, we are interested in the enhanced red- 
shift depth of background galaxies magnified by fore- 
ground clusters and averaged over large cluster samples 
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recently identified within the SDSS survey. We measure 
the mean redshift of BOSS galaxies behind SDSS clus- 
ters (lYork et al.ll2000t ISchlegel et all2009HDawson et al.l 
I2013IV using four different catalogues of optically de- 
tected clusters found from very different independen t 
methods, incl uding the maxBC G dKoester et al. l [20071. 
the GM BCG (IHao et al.ll2010h . the ISzabo etldlMrl 
and the iWen et al.l (|2012f) samples. We compare the 
depth magnification to the expected lensing signal us- 
ing mass-richness relationships calibrated from X-ray and 
weak shear measurements. In addi tion we use the MCXC 
X-ray cluster sample assembled by lPiffaretti et all ()2011l ) 
from a number of independent ROSAT observations. 

In Section §2 we present the data sets used in this 
study; the BOSS sample, the four galaxy cluster catalogs 
and the X-ray sample. In a third section, we describe our 
model, and in a fourth section we present our measure- 
ment technique, followed by our results presented and 
discussed in Section §5. In Section §6 we present a brief 
error analysis and we conclude in Section §7. Everywhere 
we assume a flat ACDM cos mology with fl„?, = 0.258, 
A = 0.742, and h wo = 0.72 (|Hinshaw et al.ll2009| ). All 
masses are expressed in unit of Mq . 

2. DATA 
2.1. The lens samples 

The Sloan Digital Sky Survey (SPSS. l^r¥eT^[200fl 
is an optical (ugriz filters) photometric and spectro- 
scopic survey which represents the largest probe of the 
local universe to date. In addition to its unprecedented 
volume observed, the success of the SDSS project is 
based upon an accurate photometric calibration, an in- 
genious synergy between photometric and spectroscopic 
observing strategies, and a handy and well documented 
database. In the published literature, several methods 
have been used to construct cluster samples from the 
SDSS database. In this study, we thus focus on the fol- 
lowing four publicly-available cluster catalogs based on 
optical identifications. 

First, the maxBCG cluster catalog (Kocstcr et al. 
2007) is a cluster sample constructed from 7 500 deg 
of photometric data in the SDSS. To detect candidate 
clusters, the method identifies a local overdensity of clus- 
ter members along the red- sequence. The richness N200 
is defined as the number of galaxy members which are 
brighter than 0.4L* and lying on the red-sequence in- 
side the radius of r2oo, within which the mean interior 
density is 200 times the critical density of the universe. 
The public maxBCG catalog consists of cluster candi- 
dates with 10 members or more. Since the maxBCG 
sample was constructed from an earlier SDSS data re- 
lease, the probed area does not fully overlap with the 
final SDSS coverage in the Galactic south cap, where 
BOSS galaxies have been targeted. Therefore, not all of 
the 3 000 deg 2 of BOSS data available to date overlap 
with the maxBCG cluster sample. Keeping all galaxy- 
cluster pairs within one deg, the subsample used in this 
study includes 6 308 clusters; 

Second, more recently, IHao et all (|2010t ) proposed a 
slightly different optical-based detection method, the 
Gaussian Mixture Brightest Cluster Galaxy (GMBCG), 
which relies on both the presence of the brightest cluster 
galaxy (BCG) and the red-sequence to identify clusters. 
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Among the differences with maxBCG cluster algorithm 
finder is the way to estimate the cluster redshift, where 
for GMBCG it is estimated solely from the BCG photo- 
metric redshift ( or spectroscopic r edshift when available) . 
It is claimed bv lHao et all (|2010f l that the GMBCG cat- 
alog is volume limited out to z = 0.4; however, to ensure 
no spurious correlation between false cluster detection 
and background sources, secure foreground-background 
sample separation is necessary. For this reason we limit 
the GMBCG sample to the range 0.1 < z < 0.3, taking 
into account the typical photometric redshift error. This 
redshift selection, in conjunction with rejection of clus- 
ters lying outside the BOSS area, leads to a total of 4 631 
clusters with richness N200 greater than 10. Although an 
improved weighted richness estimator is available, we em- 
ploy the scaled richness estimator, which is defined in the 

same way as for the maxBCG catalog. 

Third, an alternative method of iDong et al.1 (2008) 
does not rely on the red-sequence, but on the peak lo- 
cations in the likelihood map generated from the convo- 
lution of the galaxy distribution in redshift space with 
aperture matched filters. These filters are constructed 
from the assumed cluster density profile and galaxy lu- 
minosi ty functions. This method has been applied to the 
SDSS (ISzabo et aLllMtl . and the overlap with BOSS al- 
lows us to use 5 646 galaxy clusters below redshift 0.3. A 
noticeable feature of this catalog is that the richness is 
computed as the sum of the luminosity of all galaxies 
above 0.4L* divided by L+, and is therefore not directly 
comparable to the definition used by the two previous 
data samples. The publicly available version includes all 
clusters with an estimated richness A200 of 20 or higher. 
In the rest of the paper, we will refer to this sample as 
"AMF" . 

Fou rth, we use the catalog produced by iWen et aTl 
(|2012| ). In this method, a friend-of- friend algorithm is 
applied to luminous galaxies using a linking length of 
0.5 Mpc in the transverse separation and a photometric 
redshift difference within 3er along the line-of-sight di- 
rection. The center is assumed to be the position of the 
BCG, identified from a global BCG sample, as the bright- 
est galaxy physically linked to the candidate cluster. In a 
similar fashion to the AMF catalog, the richness is com- 
puted from the total luminosity of all galaxies brighter 
than 0.4L*, in units of L±. The sample is limited to a 
richness threshold of 12. The algorithm was applied to 
the latest SDSS-III data release, resulting in an increase 
of factor two in area. With a larger overlap with BOSS, 
especially in the south Galactic region, we obtain 15 112 
clusters within one degree off the edge of the background 
BOSS sample. In the following, this sample will be re- 
ferred to as "WHL12", and the richness denoted A200 
for consistency with the notation adopted for the AMF 
catalog. 

Finally, we complete our set of foreground cluster cat- 
alogs with the meta-catalo gue of X-ray detected clus- 
ters (MCXC) assembled by IPiffaretti et al.l (|2011[ ). This 
data set is a compilation of X-ray detected clusters from 
a number of publicly available data collected by the 
ROSAT satellite. The sample extends to redshift ~ 0.6 
with a wide range of masses. Here we relax the maximum 
redshift range to z = 0.35 as the majority of clusters have 
secure spectroscopic redshift measurements and do not 
suffer from false detection caused by projection effects. 



The number of clusters in common over the DR9 (see 
below) released BOSS area is currently 158. 

2.2. The Source Sample 

Our background sample is extracted from the first data 
release "DR9" of the BOSS spectroscopic survey. BOSS 
aims at measuring the scale of baryon acoustic oscilla- 
tions at redshift z = 0.5 as a sensitive cosmological test. 
The first data release, covering 3 000 deg 2 , was made 
publicly available in summer 2012 through the SDSS 
III/BOSS DR9 data releas<E 

We use the "CMASS" spectroscopic sample for which 
the targets were selected in the range 17.5 < iab < 19.9, 
using color selection techniques to ensure homogeneous 
sa mple in mass between r edshifts 0.43 and 0.7. We refer 
to iMaraston et al.l (|2012t ) for a complete study of the 
galaxy stellar mass function and the verification of the 
homogeneity of the data across the full redshift range. 

The BOSS DR9 data release comprises about 3 000 
deg 2 in the south and north galactic areas. We further se- 
lect only primary spectroscopic galaxies, and remove all 
galaxies with uncertain redshift measurements, by im- 
posing the flag zwarning=0. The total number of galax- 
ies used as background sources is 316 220. 

Fig. [T] illustrates the data coverage used in this study. 
The foreground cluster sample is shown in blue and 
the background BOSS galaxy sample in red. We show 
the maxBCG clusters as an example, and note that the 
WHL12 sample has a better overlap with BOSS in the 
south Galactic part. The two other catalogues have a 
similar coverage compared to maxBCG, and the MCXC 
has a full coverage over the BOSS area but is not spatially 
homogeneous because it is drawn from several indepen- 
dent X-ray surveys. We summarize in Table [T] cluster 
sample properties. 

3. MODEL 

In this section we outline the formalism that describes 
the systematic effect of lens magnification modifying the 
source selection function, by foreground galaxy clusters. 
Lensing magnification is caused by both isotropic and 
anisotropic focusing of light rays due to the presence of 
massive foreground objects acting as gravitational lenses. 
The former effect is described by the convergence, k(t) = 
£(r)/£ C rit) the projected mass density E(r) in units of 
the critical surface mass density, 

c 2 D s 

where Di, D s , and Di s are the proper angular diameter 
distances from the observer to the lens, the observer to 
the source, and the lens to the source. The latter effect is 
due to the gravitational shear 7(7-) = 71 + 1^2 with spin- 
2 ro tational symmetry (see iBartelmann fe SchneideH 
[200l . 

Since gravitational lensing conserves surface bright- 
ness, the apparent flux of background sources increases 
in proportion to the magnification factor. This shift in 
magnitude implies that the limiting luminosity at any 
background redshift lies effectively at a fainter limit given 
by Ln m (z) / fi(z) , hence increasing the surface density 

5 http : //www. sdss3 . org/dr9/data_access/ 
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• BOSS galaxies 

• maxBCG clusters 




Figure 1. Data coverage of BOSS galaxies (blue) and maxBCG clusters (red) in equatorial coordinates. We only display the clusters 
selected to be within one deg of BOSS area used in this study. The dashed line is the Galactic equatorial plane. 



of magnified sources behind foreground lenses. On the 
other hand, the number of background sources per unit 
area decreases due to the expansion of sky area. These 
two effects compete with each other, and the effective 
variation in the source number density n e g, known as 
the magnification bias (BTP95), depends on the steep- 
ness of the source number counts as a function of the flux 
limit F. 

For background sources at redshift z, the magnification 
bias is expressed as 



n eff [> F(z)] = 



1 



> 



F(z) 



(2) 



where no is the unlensed number density of background 
sources, L the limiting luminosity of the background 
sample, and /i the magnification, 



I 



M = 



(l-«)2-| 7 |S 



(3) 



For the case of weak gravitational lensing, the magni- 
tude shift induced by magnification is sufficiently small, 
so that the source number count can be locally approx- 
imated by a power law at the limiting luminosity. This 
simplifies Eq. [2] to 



n eS (z) = n 



0(z,L)-l 



no{z) , 



(4) 



with (5 the logarithmic slope of the luminosity function 
$ evaluated at the limiting luminosity: 



P(z,L) 



dln$[L',z] 



din 27 



(5) 



We retrieve the well-known result that if the count slope 
(3(z,L) is greater than unity, the net effect of magnifi- 
cation bias is to increase the source number density, or 
decrease otherwise (BTP95). If the count slope is unity, 
the net magnification effect on the source counts van- 
ishes. In the strict weak-lensing limit, the magnification 
bias is directly related to the projected mass distribution 



as n eS (z)/n {z) » 2[P(z,L) - 1]k(z)R 

The limiting magnitude and the count slope vary with 
redshift, so that the integrated magnification-bias effect 
will translate into an enhancement in mean source red- 
shift as 

_ / n cS (z) zdz 



J n cS (z) dz 



(6) 



The limiting luminosity can be computed from the ap- 
parent limiting magnitude for a given survey as 



- 2.51og 10 L(z) = i AB 



5.o gl0 |M_ K(i!) m 



where zab = 19.9 is the limiting magnitude of BOSS, cLl 
the luminosity distance, and K, the if-correction: 



K(z) = 2.5(1 + z) + 2.5 log 



10 



L(X e ) 

L(X ) 



(8) 



In our analysis, we restrict the source redshift range to 
0.43 < z < 0.7, where the BOSS target sele ction is estab- 
li shed to be unifo rm (|Dawson et al. 1 120 131 ). We assume 
a Schcchtci] (|1976[) luminosity function with the parame- 
ters measured in the VI MOS VLT Deep S urvey (VVDS, 
iLe Fevre et aII[2005T ) by lllbert et all (pool ). The authors 
provide a detailed measurement of the luminosity func- 
tion in all UBVRI optical bands of the survey over a 
wide redshift range, 0.2 < z < 2.0. This survey has the 
advantage of being considerably deeper than the BOSS 
survey, with redshift completeness to a fainter limiting 
magnitude and therefore describes better the underlying 
luminosity function, in terms of the slope as a function 
of redshift and magnitude required here for our lens- 
ing predictions. To minimize the uncertainties on the 
^-correction, we take the rest-frame T^-band luminosity 
function measured at redshift ^0.5, which best matches 
the redshifted galaxies observed in the i— band at redshift 

6 In the present analysis, we have used the full expression with- 
out the weak-lensing approximation. 
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0.5, so that this choice allows us to neglect the second 
term in Eq. [8] 

Finally, to account for luminosity evolution as func- 
tion of redshift, we employed the parametrization of 
iFaber et all (|2007[ ). such that our assumed Schechter pa- 
rameters are: 

M* = -22.27 - 1.23 X (z — 0.5) (9) 
a = -1.35. (10) 

Note that no precise normalization of the luminosity 
function is required for our analysis but only the gra- 
dient of the logarithmic slope /3 at the limiting luminos- 
ity as a function of background redshift, which has been 
well determined from the VVDS, and confirmed by other 
deep probes of simil ar volume such a s the COMBO-17 
photometric sur vey dWolf et ail I2003D or DEEP2 spec- 
troscopic survey (jWillmer et al.l l2006). In particular, the 
M+ parameter measured in the B— band is found to be 
in excellen t agreement amon g those three deep surveys 
at z = 0.5 (|Faber et al.ll2007l Fig. 7). 

The VVDS field-of-view is one deg 2 , comprising 11 034 
redshift measurements, allowing us to determine the re- 
spective levels of Poisson uncertainty and cosmic vari- 
ance into the Schechter parameters we require for our 
predictions. We quantify in Sec. 16.61 the impact of these 
assumptions on the size of the mean redshift depth en- 
hancement predicted. In particular we show that these 
model uncertainties arising from the parameterisation of 
the luminosity function are substantial compared to the 
detected signal, however smaller than our statistical er- 
rors on the measurements from th e current DR9 release . 

To date, a number of studies (Umetsu et al.l l2011bl : 
IHildebrandtlTal1l27)lirFord et al.ll20ll " have measured 
a magnification-bias signal by assuming an effective 
single-plane source redshift for a given background sam- 
ple and comparing their number density to a random 
sample. Here we do not need to approximate the depth 
as we have precise spectroscopic redshift measurements 
for all individual sources, providing a direct estimate for 
the enhancement of mean source redshift behind clusters 
with respect to that of the total sample. 

Our method developed here has two main advantages 
over the standard magnification-bias because we are sim- 
ply measuring an averaged depth enhancement of the 
background galaxies, rather than the lensing induced 
change in the surface density of galaxies. The latter 
effect requires careful correction for screening by fore- 
ground/cluster galaxies. Secondly, since we probe a rel- 
atively narrow redshift range of galaxies having a well- 
understood selection function than for typical faint back- 
ground source counts used in the magnification bias work 
located at z ~ 1 — 3. Therefore, our method enables in 
practice a relatively better understood determination of 
the statistical lensing effect for the foreground clusters 
we are examining. 

To quantify and characterize the cluster mass distri- 
bution, we compare the observed cluster lensing pro- 
files with analytic spherical halo densit y profiles. In 
the p resent study, we consider (1) the iNavarro et al.l 
(|1997l hereafter NFW) and (2) the singular isothermal 
sphere (SIS) models. The former is a theoretically- 
and observationally-motivated mod el of the internal 
structure of cluster-sized halos (e.g. lOkabe et al.l 120101 : 



lUmetsu et aT1l2011a|) . p(r) oc (r/r s ) _1 (l + r/r s )~ 2 with 
r s the characteristic radius at which the logarithmic den- 
sity slope is isothermal. The latter model provides a 
simple, one-parameter description of isothermal density 
profiles, p(r) oc r -2 . 

The two-parameter NFW model can be specified by 
the degree of concentration, C200 = ^200 / r s > and the halo 
mass, 



Mo 



y 200 pcrit(z)r: 



3 

200' 



(11) 



the total mass enclosed within a sphere of radius ^200, 
within which the mean interior density is 200 times 
the critical mass density at the cluster redshift, p C rit(z). 
Here we employ the mean concentration-mass relation of 
iBhattacharva et al.l (|2013| ). derived from ACDM cosmo- 
logical iV-body simulations covering a wide halo mass 
range of 2 x 10 12 -2 x 10 15 M S /iT 1 and a wide red- 
shift range of z — 0-2. We use their fitting formula 
for the full-sample relation c 2 oo(A^20o, z). The NFW 
pro file is thus specified by M 2n o alone. The use of 
the IBhattacharva et all (|2013f ) relation is motivated by 
recent detailed clu s ter lensing work (iCoe et al.l l2012t 
lUmetsu et all l2012t lOkabe et all I2013D finding a good 
agreement with their predictions for high-mass clusters 
(M 200 ~ 10 15 M Q ) at z = 0.2-0.4. We employ the radial 
de pendence of the proj ected NFW lensing profiles given 
bv I Wright fe Brainerdl (pOOOh . 
For the SIS model, the magnification is given by 



Psis(0) 



1 - Oe/0 



(12) 



with 9e = 4:ir(<7 v /c) 2 Di s /D s the Einstein radius. Here, 
a v is the one-dimensional velocity dispersion related to 
the halo mass, 



-200p crit (2)M 200 



G 3 



1/6 



(13) 



Finally, the magnification factor is plugged into Eqs.0] 
and[5l so that the depth magnification is computed as a 
function of the physical distance at the cluster redshift 

as 

z{r) 



1. 



^total 



(14) 



The total effect of magnification on the redshift distri- 
bution n(z) of BOSS source galaxies is shown in Fig. [5] 

4. MEASUREMENTS 

Given the sparse source density around foreground 
clusters, we must measure the signal around stacked clus- 
ters to increase the signal-to-noise ratio. 

To ease the computation of the mean redshift of BOSS 
galax ies around clusters, we e mployed a modified version 
of the lDavis fc Peebles! ()1983l ) two-point cross-correlation 
estimator: 

w(r)= A (15) 

where L is the lens cluster sample, S the source BOSS 
galaxy sample weighted by their redshift and R the un- 
weighted BOSS galaxy sample. Pairs are computed at 
the redshift of the cluster as function of physical scale. 
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Figure 2. Redshift distribution of BOSS background galaxies 
lensed by foreground SDSS clusters (in blue: dashed line for our 
model; dotted line with error bars for the AMF measurements) 
compared to the unlensed distribution (black solid line), and av- 
eraged within r = 0.2 Mpc from the cluster center. The model 
assumes an effective cluster mass of M200 = 1.81 X 10 14 M o at 
z = 0.2. The observed enhancement in mean source redshift is 
Sz ~ 0.01. 



We can write the previous equation explicitly: 



ni(r),n s (r) 



X Zn 



so that 



and 



LS(r)= Yl 1 

i=i,j=i 
LR(r)=n\(r) x n s (r) 



j-^{r) = z(r) 



a- x = — y 1 



X Zj — Ztotal • 



(16) 
(17) 

(18) 
(19) 



i=X,j=l 



where the subscripts I and s denote the lens and source 
samples, respectively. Equation IT51 then reduces to 



w(r) 



z ( r ) 

•Ztotal 



1 = 5 z (r) 



(20) 



as defined by Eq. [TH 

In practice we perform the measu rements using the 
software Swot (|Coupon et al.l [20121 ) . a fast two-point 
correlation code optimized for large data sets. The algo- 
rithm makes use of large scale approximations, tree-code 
structured data, and parallel computing to considerably 
accelerate pair counting. Correlating ~ 300 000 back- 
ground objects with ~ 6 000 foreground objects takes 
about five minutes on a desktop computer. 

Error bars are estimated by generating foreground 
cluster positions over the BOSS area: using 20 000 
points with random projected positions on the sky and 
a volume- weighted random redshift in the range 0.1 < 
z < 0.3 , we compute w(r) in a similar fashion as for the 
data samples, using the same binning. We repeat the 
process 100 times and compute the dispersion of the en- 
semble. Error bars for the signal of each cluster sample 



are further multiplied by y/20 000 /Clusters to account 
for Poisson error. Additionally, the significance of the 
detection for each sample is computed using the rescaled 
covariance matrix. 

For optically detected cluster samples, we compare our 
measurements to the expected theoretical signal by first 
converting the richness to the mass, using the maxBCG 
mass-richness r elationship calibra ted from X-ray and 
weak lensing bv lRozo et al.l (|2009l ): 



10 14 M, 







exp(0.95) 



200 



40 



(21) 



and we further con vert M500 into M200 using the 
method described in IHu fc Kravtsovi (I2003D. and the 
concen tration-mass relationship of iBhattacharva et al.l 
(|2013h . We employ the same relation for GMBCG, as 
the richness definition is identical to maxBCG. 

As for the WHL12 sample, the authors chose a different 
definition for the richness (see Sec. I2.ip . but also provide 
an independ ent X-ray /lensin g calibrated mass-richness 
relationship (jWen et al.ll2010( ): 



21 11 ) 



10 



-1.49 A 1.17 



(22) 



Since the richness definition assumed for AMF is similar 
to that of WHL12 we use the same parametrisation for 
these two samples. 

For the X-ray cluster sample, an estimation of the A/500 
mass is directly provided, however due to the low an- 
gular resolution of ROSAT, X-ray mass estimates may 
be underestimated. To correct for this, we matched the 
MCXC clusters with t he local cluster sample presented in 
Vikhl inm et al.l (|2009l ). and compared the value of M 50 o: 



AL 



500 



10 



-0.951 i\,fl.067 
IV1 500 ■ 



(23) 



Then we simply convert M 500 into M200 as described 
above. It is found that the MOX C-based value is about 
25% smaller than those derived bv lVikhlinin et al.l (|2009l ) 
at the effective mass of our subsamplc M200 ~ 5 x 
1O 14 M . 

Finally to perform a strict comparison between the 
stacked measured signal and the prediction, we compute 
the composite halo-mass theoretical signal: 



1 



5z{r) = -^^£ z ,i(r,M 2 oo) 

1 „- 1 



(24) 



5. RESULTS 



We show in Fig. [3] our measurements for the four 
maxBCG, GMBCG, AMF, and WHL12 SDSS cluster 
samples. Similarly, Fig.|4]shows our results for the X-ray 
MCXC cluster sample. 

As described in Sec. [U the errors are obtained from 
the standard deviation of 100 subsamples with 20 000 
randomized cluster positions over the BOSS area, and 
rescaled for each individual sample depending on the 
number of clusters to account for Poisson uncertainty. 
The detection significance, shown in Table Q] for each in- 
dividual sample, is computed for all data points with 10 
cluster-galaxy pairs or more, using the covariance matrix 
of the random samples up to r — 10.0 Mpc. The result- 
ing covariance matrix is displayed in Fig. [5] the higher 
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r (Mpc) r (Mpc) 

Figure 3. Mean redshift increase of BOSS galaxies measured around SDSS clusters from the maxBCG (top left), GMBCG (top right), 
AMF (bottom left), and WHL12 (bottom right) samples. The red dashed (composite) and black solid lines are NFW profile predictions 
assuming a mass-richness relationship calibrated from X-ray and weak lensing measurements taken from the literature. We also show the 
SIS profile prediction (dotted dashed line). Error bars were computed by calculating the background mean redshifts radially around 20,000 
randomized foreground positions over the BOSS area and normalised to the cluster sample size. 
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Figure 4. Same as Fig. [3] but for the X-ray MCXC sample. 



correlation at large scales translates into the fact that the 
same background galaxy may be used for several cluster- 
galaxy pairs. On the other hand, there is little correlation 



at smaller scales (< 1 Mpc) due to the relatively sparse 
distribution of foreground clusters (approximately, two 
per sq. deg). 

For each case, the red-dashed line represents the 
composite-NFW lensing signal expected for the respec- 
tive cluster sample, obtained assuming a mass-richness 
relationship calibrated from X-ray and weak-lensing mea- 
surements taken from the published literature: given by 
IRozo et all (120091) for maxBCG and GMBG (Eq.[H) and 
bv iWen et alTpOll ) for AMF and WHL12 (Eq. [24]), 
whereas the M200 masses for MCXC have been trans- 
lated from their M500 values assuming the NFW form. 

For comparison, we also plot the expected lensing sig- 
nal obtained assuming a single effective mass (M200), for 
the NFW (black solid line) and SIS (black-dotted dashed 
line) profiles. At small scales, there is a small difference 
between the composite versus single NFW profiles, which 
however is not significant due to the large statistical er- 
rors. As for the SIS case, overall, it is likely that the 
lensing signal at small scales is overestimated. Here, it 
is reassuring to see that the NFW model provides a bet- 
ter description of the observed lensing signals because the 
lensing-based mass-richness relationships were calibrated 
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Table 1 

Properties of cluster samples. 



Sample 


-^clusters 


Effective M 2 oo mass [xlO 14 M ] 


Redshift 


Detection method 


Significance (cr) 


Reference 


maxBCG 


6 308 


1.48 


0.1 < 2 < 0.3 


red sequence 


4.8 


1 


GMBCG 


4 631 


1.40 


0.1 < 2 < 0.3 


red sequence + BCG 


2.8 


2 


AMF 


5 646 


1.81 


0.1 < 2 < 0.3 


aperture matched filter 


4.9 


3 


WHJ12 


15112 


1.18 


0.1 < 2 < 0.3 


photo-z 


3.5 


1 


MCXC 


158 


5.00 


0.1 < 2 < 0.35 


X-ray 


2.8 


5 



References. — (1) Koostc r ct al. (20071); (2l lHao etal] 120101') ; O USzabo et~aT1 H201M : f41 IWen et~all 11201211 ; (5jPiffarctti ct ""aT1 H2011I1 . 
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Figure 5. Correlation matrix of the covariance matrix estimated 
from the 100 subsamples with randomized cluster positions in the 
BOSS area. The result shown here was computed using 20, 000 
random points for each of the 100 subsamples. 



assumin g the NFW form in all cases for the maxBCG 
relation 1 Johnston et "all 120071: iMandelbaum et all 12051 
ISheldon et alj|2009r). an d most cases for the WHL12 re- 
lation ()Wen et al.Tl2010L and references therein). 

These results demonstrate the robustness of our mea- 
surements and appear to be in excellent agreement with 
previous lensing studies. Although the large statistical 
errors prevent us from accurately testing both the ampli- 
tude and slope of the mass-richness relationship for the 
current BOSS sample, it is interesting to compare the re- 
sults between the different cluster samples. With a detec- 
tion significance at the nearly 5ct level, the maxBCG and 
AMF samples seem to give the highest-confidence detec- 
tions among our cluster samples. The AMF catalog cat- 
alogue has the highest minimum richness cut (A200 > 20) 
but the averaged signal appears to be relatively higher so 
that the lower number of clusters above this richness cut 
is almost compensated by a greater lensing signal due to 
their higher cluster masses. 



On the other hand, the measured lensing signal in the 
WHL12 sample is shown to be relatively low compared 
to the model prediction based on the same mass-richness 
relationship, despite its larger sample size and greater 
overlap with current BOSS sky area. This could be a 
consequence of a relatively higher impurity of the cluster 
sample, due perhaps to a higher rate of false detections, 
which will lead to an underestimation of the geometric 
lensing signal as found here. Similarly, we obtained for 
GMBCG a relatively low level of detection s i gnific ance, 
2.4(7. It is interesting to note that lWen et al.l (|2012f ) find 
a lower matching rate of their WHL12 clusters with GM- 
BCG, compared to that with maxBCG. This could indi- 
cate that these two samples have a higher level of con- 
tamination by false cluster detections. Another inter- 
esting point is that the cluster redshift scatter between 
WHL12 and GMBCG can be as high as several times 
the phot ometric redshift errors, as shown by their com- 
parison (|Wen et al .112011 Table 8). A more detailed and 
careful analysis would be certainly required to confirm 
such a conclusion. We shall come back to this issue and 
its consequence for magnification bias contamination in 
Sec. EU 

The MCXC X-ray sample has the largest statistical er- 
rors due to the small number, 158, of clusters (Fig. [4]). 
However, it is worth noting that there is a likely excess 
of the lensing signal at r ~ 2-10 Mpc (> r2oo) with re- 
spect to the NFW predictions, which cannot be explained 
by simply increasing the halo masses. In the context of 
ACDM, this large-scale excess signal can be naturally ex- 
plained by the two-halo term contribution due to large- 
scale structure associat ed with the central clusters (see 
IQguri fc Hamanal 120111 ). We note that, in contrast to 
the gravitational shear, magnification is sensitive to the 
sheet-like mass distribution, and therefore can be used as 
a powerful tool to probe the two-halo term in projection 
space. 

6. ERROR ANALYSIS 

6.1. Sources Sample Variance 

As described in Sec. 2J our error estimate simply re- 
flects the statistical variation of the mean redshift of 
sources across the BOSS area. This estimator primarily 
provides a measure of the source sample variance given 
a sample of foreground clusters, because it measures the 
fluctuations of the number counts across the field. We 
display the magnification signal as measured for AMF 
(for which we find the highest level of detection signifi- 
cance) in the top panel of Fig. [SJ compared to the mean 
signal from the random samples shown in the bottom 
panel. The random-sample signal is consistent with null 
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Figure 6. Measured signal for the AMF sample (top) compared 
to the mean signal around random positions drawn in the BOSS 
area (bottom). The error bars show the dispersion for 100 samples 
with randomised cluster positions, mainly accounting for the source 
sample variance around the cluster sample. In blue we display the 
expected NFW signal as in Fig|3] where the shaded area represents 
an estimate on our model uncertainties due to the assumptions on 
the luminosity function. 



detection, indicating that the 100 random samples are 
sufficiently large for realistic estimates of the errors. 

6.2. Spectroscopic Redshift Errors 

The errors on individual source redshifts might also af- 
fect our statistical lensing measurements if only a small 
number of cluster-galaxy pairs are used at small scales 
( ;$ 0.05 Mpc). However, the mean spectroscopic error 
is Az = 1.3 x 10~ 4 , so that our measured signal is tow 
orders of magnitude larger than that. Hence, we con- 
clude that this should not have a significant impact on 
our measurements. 

6.3. Photometric Calibration 

The primary requirement of our analysis method is an 
homogeneous target selection across the survey field. In 
particular, a spatial variation of the photometric cali- 
bration would affect both the cluster detections and the 
BOSS target selection. In our analysis, we can rely on 
the a ccurate and homogeneous calibration o f the SDSS 
data (|Fukugita et al.lll996t iGunn et al.lll998D . We have 
checked that a variation of a few 0.01 magnitude would 
cause a mean redshift change less than Sz = 5 x 10~ 4 , 
which is again significantly smaller than the measured 
signal. 

6.4. False Cluster Detections 

Perhaps the most important source of systematics in 
magnification bias methods originates from the poten- 
tial contamination by physically associated lens-source 



pairs that create a spurious correlation signal indistin- 
guishable from the lensing signal. Since we use spec- 
troscopic redshift information for the background, low- 
redshift contamination by background sources can be ne- 
glected here. Hence, the only possible source for such 
systematics would arise from erroneously-detected dis- 
tant clusters. Let us consider a galaxy overdensity (such 
as rich clusters or filamentary structure) located in the 
source range 0.43 < z < 0.7 whose member galaxies were 
misattributed to foreground galaxies with lower photo- 
metric redshifts, which can be incorrectly detected as a 
low redshift cluster. Since cluster-sized overdensities can 
have a significant impact on the redshift distribution, the 
resulting mean redshift could be noticeably different from 
that of the whole sample. 

The maxBCG and GMBCG samples would have the 
lowest chance of being affected by this effect, because 
their detection algorithms rely on bright red galaxies, for 
which wc expect a very low rate of photometric redshift 
outliers (defined as being off by a few sigmas) due to a 
clear measurable Balmer break in the redshift range con- 
sidered here, and the high rate of spectroscopic redshifts 
for these objects. On the other hand, for the WHL12 and 
AMF samples which make use of the full galaxy popula- 
tion with photometric redshifts, such misidentification is 
more likely to happen. 

Interestingly, Fig. 4 in lOvaizu et al.l ([2008) shows a 
photometric-to-spectroscopic redshift comparison of the 
CC2 sample that was used for AMF. A closer look at the 
published figure suggests that very few galaxies (one can 
count only six galaxies out of ~ 10 000 - less than 0.1%) 
scatter from the redshift range [0.43 : 0.7] to [0.1 : 0.3]. 
Taking into account the increased scatter in overdense 
regions, a threshold richness of ~ 10 in cluster catalogs 
implies that only those clusters which are richer than a 
few hundreds would be large enough to cause such a false 
detection at low redshift. 

Also rare, but more likely to occur than in the previous 
case, is massive filaments aligned with the line-of-sight, 
which can lead to enhanced low-z contamination. This 
might cause a false cluster detection and also impact 
the mean redshift esti mation of the backgr ound sample. 
Again from Fig. 4 in lOvaizu et al.l (|2008T ) . it is clearly 
seen that the scattered galaxies lie below the mean red- 
shift of BOSS galaxies. This means that the expected 
direction of this effect would be to decrease the mean 
redshift and bias low our measurements. This could con- 
stitute a plausible explanation for the low signal level 
measured from WHL12, since they employ photometric 
redshifts in their detection algorithm. 

It is worth noting that such large overdensities would 
be easy to spot and mask out in a spectroscopic back- 
ground sample, when the statistical significance of the 
BOSS survey improves. 

6.5. Cluster Miscentering 

Misidentification of cluster centers is another poten- 
tial source of systematic errors for cluster lensin g mea - 
surements at small scales. Recently. iGeorge et al.l ()2012f ) 
studied the impact of the choice for the cluster center on 
the measured lensing signal based on X-ray galaxy groups 
detected in the COSMOS field. Their findings show that 
at small scales (< 0.1 Mpc), choosing the BCG over the 
most massive galaxy close to the X-ray peak position 
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(their center definition yielding the best result) reduces 
t he signal by abo u t 20% within 75 Kpc. 

Uohnston et al.l (|2007f) demonstrated that the lensing 
convergence k, which is locally related to the magnifi- 
cation to first order, is less affected by cluster miscen- 
tering than the weak shear, and that smoothing due to 
miscentering effects nearly vanishes at twice the typical 
positional offset from the cluster mass centroid. This 
indicates that our results would not be affected by the 
miscentering effects beyond a radius of r = 0.15 Mpc. 
In our present analysis, the statistical uncertainty is too 
large to be able to estimate the degree of miscentering 
or the likely level of correction from our data. Still, it 
is worth noting that our measurements show the lens- 
ing signal decreases at small scales for three out of the 
four samples, although these cluster samples are not en- 
tirely indep endent with sizea ble proportions of clusters 
m common (|Wen et al.ll20l"l . 

6.6. Model Uncertainties 

In this work, while we attempt to demonstrate the 
feasibility and great potential of this new method, we 
are limited by statistical uncertainties and thus unable 
to constrain both the amplitude and slope of the mass- 
richness relationship. When BOSS is completed, it will 
increase the respective size of background sample and 
the large area of sky covered - which means we can also 
increase the size of the foreground cluster samples by a 
factor of three - leading to an anticipated factor-three 
increase in the signal-to-noise ratio. This will then al- 
low us to subdivide cluster samples into bins of rich- 
ness, thus complementing the standard shear-based mass 
measurements. A further substantial improvement is ex- 
pected in near future with bigBOSS and the NIR spectro- 
scopic survey of EUCLID, which will increase the sam- 
ple size of background galaxies by one order of magni- 
tude. More importantly, EUCLID will obtain a suffi- 
ciently large number of spectra of background sources at 
z > 1, pushing the foreground cluster sample to higher 
redshift than shear measurement will allow. Similarly, 
the ambitious Su baru Prime Focus Spectrograph (PFS, 
ISugai et al.l I20T1) survey with its powerful NIR capa- 
bility promises to provide an unprecedented large high- 
redshift sample of spectroscopic redshifts. 

To fully exploit the greater statistical power in future 
surveys, it is critical to control systematics. In particular, 
the luminosity function slopes of source samples must be 
precisely known, in order to accurately extract the lens- 
ing signal from observations. To do this, a homogeneous 
selection of background sources is crucial. 

To compute theoretical predictions, we have implicitly 
assumed that mass is proportional to luminosity when 
constructing the above composite cluster lensing mass 
profiles. In details we must expect that the mass-to-light 
ratio vari es as a function of m ass, redshift, and galaxy 
type. In iCoupon et al.l (|2012f ). we fit the mass-to-light 
ratio as function of mass and redshift for blue and red 
galaxies, using the 30-band COSMOS data (jllbert et al.l 
I2009T ): in the mass bin 10 10 5 - 10 115 M Q , we find the 
slope of the luminosity function changes from 0.44 to 0.69 
depending of the galaxy type. This implies that, at most, 
the change in luminosity keeping a constant mass would 
be about 0.1 magnitude from redshift 0.4 to 0.7. Given 
the parametrization of the luminosity function employed 



here we can test how the slope of the number counts 
would be affected by simply shifting the limiting mag- 
nitude by 0.1. In practice, since we are concerned only 
with the relatively bright end of the galaxy distribution, 
it is equivalent to a shift in M+ parameter. 

Furthermore, two additional sources of systematic er- 
rors could have affecte d the determ i nation of the source 
luminosity function by lllbert et al.l (|2005l ): the Poisson 
err or and the cosmic y ariance. The former was estimated 
bv lllbert et ail (|2005l) , and given as Ail/* ~ 0.25, which 
is correlated with the faint-end slope a, so that this esti- 
mate may be rather pessimistic. For the later, the effect 
of cosmic variance on the count slope is more difficult to 
assess. To obtain a crude estimate of this effect we con- 
sider the value of AM t = 0.02, which was estimated by 
ILovedav et all (pOllh for the GAMA survey from jack- 
knife resampling of nine subregions each with 16deg 2 . 
Although GAMA is a lower-redshift survey compared to 
VVDS, the volume of each subregion is larger than the 
volume probed over the 1 deg 2 field-of-view of VVDS at 
z = 0.5, so that we have adopted a conservative estimate 
of 0.1. 

Finally we add in quadrature all three sources of sys- 
tematic uncertainty, and reexamine our model predic- 
tions in the two extreme cases, M* ± 0.29. In Fig. [SJ 
the resulting error region is indicated by the blue-shaded 
area which is systematically less than the statistical error 
for the current dataset. 

7. CONCLUSIONS 

Using over 300 000 BOSS galaxies we have measured 
the mean redshift of background galaxies behind large 
samples of SDSS galaxy clusters. Our results show a net 
increase of the mean redshift behind the clusters com- 
pared to that of the total sample, in line with reasonable 
expectations for the effect of lens magnification. We have 
tested four different cluster catalogs, the maxBCG, GM- 
BCG, AMF and WHL12 samples, with detection signif- 
icance ranging from 2.8ers (GMBCG) to 4.9ers (AMF), 
where the level of systematic errors is expected to be 
negligible compared to our statistical errors. In order to 
speed-up and ease the measurements of the mean redshift 
as function of physical scale around clusters, we employed 
a sophisticated code (Swot) that performs rapid paral- 
lel calculations and accurate error estimations, making 
practical the handling of future samples containing sev- 
eral millions of objects. 

Based on precise measurements of luminosity functions 
from a number of deep spectroscopic surveys, and on 
an accurate modelling of dark matter halo profiles, we 
have compared our results to theoretical predictions. We 
constructed the expected signal by summing the indi- 
vidual contribution for each cluster given its richness 
from two independently established mass-richness rela- 
tionships corresponding to two kinds of richness defini- 
tion, and each calibrated with X-ray and weak lensing 
data. For three samples out of four, the agreement with 
the NFW profile is excellent. The mean masses of the 
clusters derived from our basic stacking analysis vary be- 
tween the four large SDSS cluster samples in the range 
1.18-1.81 14 M Q , after allowing for the mass-richness re- 
lations appropriate for each cluster sample and with a 
mean radial profile consistent with the observed radial 
trend towards higher mean redshift at smaller cluster ra- 
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dius. Only WHL12 showed a marginal deviation from 
the expected level, suggesting an underestimated signal. 
Given the high density of candidate clusters in this cat- 
alog compared to the others, this effect may indicate a 
somewhat larger contamination of spurious clusters that 
dilute the lensing signal in proportion to the fraction of 
false detections. 

Further investigation of possible sources of systematic 
error will be feasible with a more thorough understand- 
ing of the source population of the BOSS survey, in terms 
of colour and magnitude selection and Eddington bias af- 
fecting the very bright end of the BOSS luminosity func- 
tion, arising from photometric error. With the current 
sample, given our relatively large statistical errors we 
conclude that the systematic errors in the measurements 
should be negligible at present but that with the comple- 
tion of the BOSS survey we will reach an increased level 
of precision that may warrant closer scrutiny of the de- 
tails of the source selection function and possible redshift 
uncertainties for the foreground cluster populations. We 
believe that the method is very robust against system- 
atic errors, even when we increase the statistical power 
in future experiments such as bigBOSS, PFS and EU- 
CLID (which will observe background sources 25 times 
as dense and over 5 times the area of this study) with 
millions of spectroscopic/grism redshifts, at which point 
this method will be powerful in its own right for defining 
cluster mass-concentration relations and the evolution of 
the cluster mass function. These results may be com- 
pared with independent estimates of magnification bias 
from faint number counts and with cosmic shear mea- 
surements for which the sources of systematic error are 
very different. 
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